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Abstract 

If there are large extra dimensions and the fundamental Planck scale is at the 
TeV scale, then the question arises of whether ultra-high energy cosmic rays might 
probe them. We study the neutrino-nucleon cross section in these models. The elastic 
forward scattering is analyzed in some detail, hoping to clarify earlier discussions. We 
also estimate the black hole production rate. We study energy loss from graviton 
mediated interactions and conclude that they can not explain the cosmic ray events 
above the GZK energy limit. However, these interactions could start horizontal air 
showers with characteristic profile and at a rate higher than in the Standard Model. 
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1 Introduction 



In this paper we explore the possibihty that the primary particles for ultrahigh-energy cosmic 
rays are neutrini interacting gravitationally with atmospheric nucleons. An obvious objection 
to this idea is that the gravitational interaction is too weak to produce any sizable cross 
section for this process. However, this point needs to be fully reconsidered in theories where 
the fundamental scale is around the TeV, as postulated in models involving large extra 
dimensions and a four-dimensional brane- world [0, 0. In these scenarios not only does the 
cross section for elastic gravitational scattering increase at cosmic ray energies, but there is 
also the possibility that the collision results into the formation of microscopic black holes. 
Both effects can dramatically increase the cross section for scattering between neutrini and 
atmospheric nucleons, and hence they may play a role in explaining the most energetic cosmic 
ray events. 

In the recent past, this possibility has been entertained in a number of papers, with 
differing conclusions ^ There has been a controversy as to the right way to perform 
the calculations, and how to implement unitarity at high energies. We hope to shed some light 
on these issues, and show that, actually, the situation is quite simple once the appropriate 
point of view is taken 0. Besides, the possibility of producing black holes in cosmic ray 
collisions needs to be addressed in detail. Once we have, hopefully, settled the terms for the 
analysis, we will turn to the actual discussion of whether the first signatures from low-scale 
unification and large extra dimensions might come from the study of ultrahigh-energy cosmic 
rays (UHECR). 

While this work was in progress, detailed analyses of the possibility of black holes forming 
at the LHC have appeared 0, H, 0. During the final stage of this work, a paper |]lOl has 
appeared which studies the production of black holes in cosmic rays, and also investigates 



their detection in horizontal air showers. Our work complements that of ref. |T0[: the latter 
focuses on the phenomenology of the detection of black holes, whereas we address in more 
detail the theoretical aspects of neutrino-nucleon scattering in TeV-gravity theories. A paper 
discussing further aspects of the detection of these showers has appeared when this work was 
ready for submission ITlI . 



2 Ultra-high energy scattering on the brane 

An essential feature of the gravitational interaction is that at center-of-mass (c.o.m.) energies 
•y/s well above the fundamental scale, the coupling to gravitational coupling grows so large 
that graviton exchange dominates over all other interactions. This is actually the case for 
atmospheric nucleons being hit by neutrini of energy ~ 10^^ GeV {^/s ~ 10^ GeV) if the 
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fundamental scale is around 1 TeV. In particular, if the impact parameter b is sufficiently 
smaller than the radius of compactification, the extra-dimensions can be treated as non- 
compact. In this regime one would be probing the extra dimensions purely by means of the 
gravitational interactions. 

Another consequence of ultra-high energies in gravitational scattering is that, to leading 
order, its description involves only classical gravitational dynamics. In particular, this means 
that we do not need any detailed knowledge of quantum gravity to perform the calculations: 
any theory that has General Relativity as its classical limit should yield the same result^. 
One can distinguish different regimes in the scattering (and we will do so below), but perhaps 
the most spectacular effect at such energies is the expected formation of black holes, via 
classical collapse, when the impact parameter is of the order of the horizon radius of the 
(higher dimensional) black hole |l5[1p|. 



n + 2 ) 



This implies that any dynamics at 6 < Rs is completely shrouded by the appearance of 
trapped surfaces: Ultrashort distances are directly probed only for energies around the 
fundamental energy scale. 

In the following we will assume for simplicity that no scale for new physics arises before 
reaching the scale for the fundamental energy M^. In particular, we assume that scales such 
as the string tension, or the tension and thickness of the brane, do not appear before that 
scale. This prevents the possibility of additional effects arising at impact parameters larger 
than the ones that give rise to black hole formation. If this is the case, then the picture for 
ultrahigh-energy scattering that we describe here should be largely universal. Nevertheless, 
stringy effects below the regime where General Relativity can be trusted may be readily 
accommodated and should not introduce large changes in our results. 

Ultra-high energy scattering in the Randall- Sundrum model has been addressed in Q, 
and the different regimes for the scattering in the present case are qualitatively the same 
as described there. At large impact parameters one does not expect formation of black 
holes, but in this case, the leading contribution to the scattering amplitude is exactly (non- 
perturbatively) calculable within an eikonal approach [l^, |13|, This is known to work 



particularly well for high energy gravitational scattering at large impact parameter |T3|, |l^ P| 



The eikonal resummation of ladder and crossed-ladder diagrams is achieved by computing 



^This has been noted often earlier, e.g., in ||T^, [T^, 14, |[ ^ ||. 
2 We define Md as in [0. 

^Loops involving only momenta of internal gravitons are suppressed by factors of l/{Mi:)b)^. 
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the scattering amplitude as 

Mis,t) = y 1 d^foe^'i-b (e 
ins 



dbbJoibq) (e^'^^^''') - l) . (2) 



This amphtude is well defined for any values of the exchanged momentum q = \f—i (t < 
since the scattering is elastic). The eikonal phase xi^s^l)) is obtained from the Fourier- 
transform to impact parameter space of the Born amplitude. Alternatively, it can be obtained 
from the defiection of a particle at rest when crossing the gravitational shockwave created 
by a second particle |12[ . 

Note that the transforms in impact parameter space are two-dimensional, since the par- 
ticles scatter in three spatial dimensions. Nevertheless, the exchanged gravitons propagate 
in the 4 + n dimensional space. Moreover, we are working at a scale where the spectrum 
of Kaluza-Klein modes is essentially continuous. In this case the Born amplitude comes out 
easily as [|l^] 

iMsorn = ^7r"/^r (l - ^) ^("t - ^e)^-^ . (3) 



1 ^2 

Xis, h) = — I —±-e^^-\MBorn , (4) 



Hence the eikonal phase 



2isJ (27r)2 

which is finite for 6 7^ 0, is xi^-, b) = (bc/b)"", where we have defined 

b- - M!:ir (-] (5) 

Having the phase x{s, b) is sufficient for numerical evaluation of the eikonalized amplitude 
d^). The result in eq. |^ can be written in terms of Meijer functions. However, it is easy to 
get simple analytical expressions for the amplitude in both regimes of qbc ^ 1 and qb^ <C 1.. 
When q ^ b~^ the phase x{S) b) yields a sharp peak for the eikonal amplitude in (0), which 
allows for an evaluation near the saddle point bg = bc{qbc/n)~^^^'^~^^^ -C b^. 



M 



saddle 



(47r)t-ir 




n+2 



n + l 



n+2 
n + l 



*This corresponds to the linearized approximation to 
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The phase (p = {n + l){bc/bs)"' is real for t < 0. Observe that the amphtude is non- 
perturbative in the gravitational coupling l/M]^^. In the limit q ^ one gets instead 

M{q = 0) =2m sblT (l--)e^ , (7) 



which is finite for n > 2. For n = 2 the real part of Ai has a logarithmic singularity 

M -Ansbl Hqb,). (8) 

Notice that also at small q the amplitude is non-analytic in the gravitational coupling. Indeed 
the amphtude at g — > is effectively described by the (Born) operator T defined in ref. 
but with an effective UV cut-off ~ b~^ on the mass of the exchanged KK modes. This cut- 
off originates from the interference with the multigraviton exchange diagrams in the eikonal 
series. 

The eikonal amplitude will be used in the next section to compute the differential cross 
section for neutrino-nucleon scattering. At the partonic level we have 

dq'^ IGvrs^ ^ ^ ^ ^ 

We can also derive the total elastic cross section from the optical theorem: 

(Tel = = 2Txbi r 1 - - COS - , (10) 

s \ nj n 

i.e., it is essentially given by the area of a disk of radius ~ fee- 
Observe that 

aeZ-S^/". (11) 

This growth of the cross section at high energy is slower than the perturbative result a ~ s^, 
and also slower (for n > 2) than the linear dependence a ~ s postulated (apparently for 
all n) in 0]. Unitarity in impact parameter space is manifest in the eikonal amplitude (0)|. 
For large impact parameter this implies as well unitarity for high partial waves. Partial 
wave unitarity at shorter impact parameter is a harder problem, and indeed, corrections to 
the eikonal amplitude are expected to become crucial. As t grows, graviton self-interactions, 
which carry factors of t associated to the vertices, increase the attraction among the scattered 
particles, and it is expected that, eventually, gravitational collapse to a black hole will take 
place. Hence the initial state is expected to be completely absorbed, but in such a way that 



^Notice that in order to achieve unitarity we have needed to perform an all-order loop resummation. As 
argued in JT^ , this is essential when considering energies above M^. This point is missed in some of the 
earlier work, such as the last reference in Note as well that the Froissart bound generalized to 
higher dimensions in does not apply since the exchanged particle is massless. 
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any short distance effects will be screened by the appearance of a horizon. Indeed as shown 
in ref.[^ the effects of the non-linearity of gravity are suppressed by a power of Rs/b, so our 
eikonal approximation should be valid for b ^ Rs and its breakdown be associated to the 
formation of black-holes. This relation between eikonal breakdown and black-hole formation 
can also be establisehd as follows. In the region b <^ be, there is a one to one correspondence 
between the transferred momentum q and the saddle point impact parameter bs. The case 
q ~ ^/s, where the (small angle) eikonal approximation breaks down, corresponds precisely 
to bg ~ Rs- Notice in passing that we can also write eq. (§) as da = 2nbsdbs, as expected 
for a classical trajectory with impact parameter bg. 

At present, the cross section for black hole production can only be estimated as the 
geometric cross section, 

~ TiRl . (12) 

with Rs as in (|T]). In this case ~ again slower than linear. 

Clearly this result cannot be very accurate. Radiation is expected to be emitted during 
the collapse, and the amount of energy that is expected to be radiated in the process can 
be a sizable fraction of the total energy (perhaps around 15 — 30%, from four- dimensional 
estimates pH), but at large enough energies it will not be able to prevent the collapse. 
This effect will tend to reduce the above value for the cross section. However, there are also 
factors which increase it, such as the fact that a black hole acts as a somewhat larger scatterer 
(40 — 75% larger radius ||2^). It seems reasonable to expect that the above expression is 
not off by any large factors. Finally, note that these black holes form through classical 
collapse. In the semiclassical instanton contribution to the nucleation of black holes 
was considered. Being a tunneling process, it is exponentially suppressed. Hence, it can be 
neglected relative to the (real time) classical collapse we are considering. 

3 Neutrino-nucleon scattering and black hole produc- 
tion 

These results can now be readily applied to neutrino-nucleon scattering at ultra-high energies. 
At impact parameters 6 < 1 GeV"^ the neutrino interacts essentially with the partons, and 
if b > Rs the eikonal approximation gives a good description of the scattering. At smaller 
distances, trapped surfaces are expected to form and the neutrino and the parton will collapse 
to form a black hole. 

In order to numerically evaluate the amplitude (0), we proceed as follows. First, we write 
it as 

iM = Aixsbl f xdxJ{xqbe){e'^'^" - 1) = AnsblMiqb,) . (13) 
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At large values of qbc we know this is well described by the simple result (^) . It is convenient 
to extract this behavior, and write the squared amplitude |7Vlp as 

2 

^ / \ _I1±1 n n + l 

\M\' = {l + (qKr) "^^ ^H<lbc). (14) 

The prefactors have been chosen in such a way that for qbc —>■ oo the function F goes to 1. 
Apart for the case n = 2 where it has a mild logarithmic singularity at qbc —>■ (see eq. (P)), 
F is 0(1) over the full range of qbc- 

For our applications it is useful to study the cross section as a funtion of the fraction y 
of energy transferred to the nucleon: 

E„ - El q^ 

y = — ^ — = — • (15) 

Fi, xs 

where x is the fraction of proton momentum carried by the parton. Summing over partons 
we have 

Here fi{x,fi) are the parton distribution functions (PDFs) (we use the CTEQ5 set extended 
to x < 10~^ with the methods in 0]). Notice that quarks and gluons interact in the same 
way. The scale /i should be chosen in order to minimize the higher order QCD corrections 
to our process. A simple, but naive, choice would he fi = q. However 1/q does not really 
represent the typical time or length scale of the interaction. As we have seen, in the sta- 
tionary phase regime, the neutrino is truly probing a distance bs ^ 1/q from the parton. 
Heuristically: the total exchanged momentum can be large, but through the exchange of 
many soft gravitons. So we believe that a better normalization is to take fi = when 
q > b^^ and /i = g if g < b^^ . The latter choice is effectively equivalent to choosing /i = 6^^ 
as at small q the eikonal corresponds to a pointlike interaction. Our choice of // is consistent 
with the fact that gravity at ultra-Planckian energies is dominated by long distance classical 
physics. Choosing fi = q would also make little sense, q can be as big as ~ ^/s ^ Md, but 
the evolution of the PDF's at > M|, cannot be simply performed withing QCD, as truly 
quantum gravitational effects (string theory) would come into play. Instead as ^/s grows 
above Md, and t/s is kept fixed but small, the impact parameter bs grows and we are less 
sensitive to short distance physics. As a matter of fact, for large enough s the total ct^n 
will be bigger than the proton area ~ (GeV)^: at higher energies the parton picture breaks 
down, the proton interacts gravitationally as a pointlike particle, and the neutrino scatters 
elastically on it. 

A useful quantity to study is the cross section integrated for y > y^. In Fig. 1 and Fig. 2 
we plot this quantity for Mo = 1 TeV and = 5 TeV, respectively. We include the cases 
with n = (2, 3, 6) and F^ = (lO^o, lO^^, lO^^) GeV. 
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Figure 1: Elastic cross section vs. minimum fraction of energy lost by the neutrino for 
Md = 1 TeV and n = 2, 3, 6 large extra dimensions. Solid, long-dashed and short-dashed 
hues correspond respectively to = 10^^ GeV, = 10^^ GeV and Ei, = 10^° GeV. 



Finally, to estimate the total cross section to produce a black hole in a neutrino-nucleon 
scattering we compute 

^= A (EM^^f^)] (17) 

where Ks is given in Eq. (|ip and \i = . A gain for the choice of scale in the PDF's the 
previous discussion applies: the Schwarzschild radius rather than the black-hole mass sets 
the time scale of gravitational collapse. Notice that in a more standard case of, say, neutrino- 
quark fusion into an elementary lepto-quark the right choice would be /i of the order of the 
lepto-quark mass. The crucial difference is that the black-hole is not an elementary object: 
its physical size is much bigger than its Compton wavelength. 

We plot in Fig. 3 this cross section versus the energy of the incoming neutrino for 
n = (2, 3, 6) and M = (1, 5) TeV. We include plots with xs > Ml (sohd) and xs > {lOMof 
(dots). These correspond to the cross sections for producing black holes with a mass larger 
than Md or IOMd, respectively. The sizeable difference between the two choices of a mini- 
mum X, indicates that the production of light (small) black holes dominates: the fast decrease 
with X of the p.d.f.'s wins over the growth oc x^''*-""'"^-' of the partonic BH cross section. No- 
tice, on the other hand, that the total elastic cross section is less dependent on the small x 
region and is dominated by x ~ 1 for the case n < 3. This is because of the faster growth 
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Figure 2: As in Fig. 1 but for Md = 5 TeV. 
oc x^/" of the partonic elastic cross section. 

4 Discussion 

We are now ready to discuss the implications of our results on the phenomenology of ultra- 
high energy cosmic rays. The first question is whether neutrino nucleon scattering at super- 
Planckian energies can explain the observed cosmic ray events with energy E > Egzk = 

5 X 10^° GeV. It is known since long ago that cosmic protons with energy above Egzk 
are damped by inelastic scattering with the microwave background photons. The relevant 
reaction isp + 7 p + it, and Eqzk is the threshold proton energy given the photon 
temperature. Because of this reaction, ultra energetic cosmic protons are brought down 
to £■ ~ Egzk within a few Mpc. Since there are good reasons to believe that the cosmic 
protons have extra-galactic origin, we should observe a sharp drop in the observed event 
rate aX E > Egzk- However, various experiments do not observe this drop at all. There 
have been several suggestions to explain that. One idea is that the primary particles for the 
UHECR are neutrini |]25|, ^, as these particles interact negligibly with the microwave 
background and are essentially undamped. However, any of these suggestions has to face the 
fact that, within the Standard Model (SM), the neutrini interact too weakly also with the 
nucleons in the atmosphere. In order to explain the ultra-GZK events by cosmic neutrini 
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Figure 3: Cross section for black hole production as a function of E^,, for Md = 1,5 TeV and 
n = 2, 3, 6. Solid and dotted lines correspond to xs > Mf^ and xs > (IOMd)^ respectively. 



one needs new physics enhancing their cross section with nucleons at high energy. In ref. 
it was suggested that, in models with TeV scale gravity, the eikonalized cross section could 
be of the right order of magnitude. However ref. [^] did not investigate the rate of energy 
loss in the eikonalized process, and, in particular, did not pay attention to its "softness". 
The production of black holes was also neglected in ref. ||^ . 

As a matter of fact, in order to determine the signal it is important to establish quan- 
titatively which is the process that dominates energy loss - whether elastic gravitational 
scattering or black hole production. It turns out that energy loss is mostly determined by 
black hole production and by scattering at ?/ ~ 1. (As we already pointed out, and as can be 
seen from comparing Figs. 1,2 with Fig. 3, the gravitational cross section at ?/ ~ 1 becomes 
comparable to a^h, though its precise value is not calculable within our linearized gravity 
approximation.) To see this, consider a neutrino travelling through a medium of density p. 
The mean free path for black hole production, at which all energy is lost to the shower, is 
Lbh = {o'bhp)~^- While travelling through the medium the neutrino also loses energy through 
the softer, but more frequent, eikonalized scattering. After travelling a distance Lbh, the 
energy fraction lost to soft scatterings with y < i/q is controlled by the quantity 

fvo da 1 pyo da 

viVo) = y-rpLbhdy= — / y—dy. (18) 
Jo dy abh Jo dy 

When 7] is less than 1 the soft scatterings play a negligible role in the transfer of energy to 
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the atmosphere. In Fig. 4 we plot t] for several cases: they all show that black hole formation 
and scattering at large y dominate energy loss. Notice that, by the discussion at the end of 
the previous section, energy loss is thus dominated by parton scatterings with ^/xs ~ Md, 
i.e. in the Planckian regime. 




1. xlO"^ 0-0001 0.01 1 
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Figure 4: Fraction 77, defined in eq. (|T8|), of neutrino energy lost to soft scatterings. Solid, 
long-dashed, and short-dashed lines correspond to E^, = 10^^, 10^^ and 10^° GeV, respec- 
tively. 

The observed showers above the GZK cut-off are all consistent with an incoming particle 
that loses all its energy to the shower already in the high atmosphere. From the above 
discussion, low scale gravity could explain these events if the mean free path Lbh for black hole 
production were somewhat smaller than the vertical depth of the atmosphere. In standard 
units, the vertical depth is measured as the number of nucleons per unit area x„ = 
1033 X A^^/cm~^ = mb~^ (where Na is the Avogadro number), so the requirement is abh > 
Xv~^ =mb. From Fig. 3 one can see that, at the relevant energies, the black hole cross section, 
however large, falls short of this requirement. In order to satisfy a^h >mb the gravity scale 
Md should be well below a TeV, which would contradict collider limits. 

Hence, we conclude that neutrino-nucleon interactions in TeV-gravity models are not 
sufficient to explain the showers above the GZK limit. Also, at present, cosmic rays do not 
appear to place any significant bounds on such scenarios. 

Nevertheless, neutrino-nucleon cross sections ai^N in the range 10~^ mb to 1 mb, like 
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in our scenario, can still lead to interesting new phenomena in cosmic ray physics, which 
may be observed in upcoming experiments. Cosmic primaries with cross section below 1 
mb can travel deep into the atmosphere before starting a shower. In particular they can 
cross the atmosphere at a large zenith angle and start characteristic horizontal air showers^. 
The horizontal depth of the atmosphere is Xh is about 36 times the vertical one, so that 
for a^N ^ .1 mb a neutrino can travel horizontally down to the interaction point. In the 
Standard Model the charged-current cross section is a^i^ ~ 10~^{E /lO^^GeV)^'^^^ mb. No 
horizontal air shower has been detected so far. However, conservative estimates of the flux of 



ultra energetic cosmic neutrini |]26| suggest that the next generation of experiments should 
be barely sensitive to neutrino cross sections of the order of the SM one. In our scenario a^N 
can be considerably bigger, so there is the interesting possibility that gravitational scattering 
and black hole production will lead to a sizeable event rate, higher than in the SM. 

The shape of the shower is probably one of the better ways to characterize these processes. 
In the SM charged-current process, a significant fraction of the neutrino energy is released 
to just one or a few hadrons from the brakdown of the target proton. The shower then 
builds up from the cascading hadronic interactions of these few hadrons. In the scenarios 
we are considering, the production of a black hole of mass M^^ ~ y/s = ^j2M^mp is followed 
by its very quick evaporation by emission on the brane ||22[ of a number of particles of the 
order of Mbh/Tbh ~ (v^/Mfl)('^+2)/("+i). For the energies we are considering this number 
can be bigger than 100. Then the shower builds more quickly than for SM processes. It 
is reasonable to expect that the shapes will differ, very much in the way that a shower 
formed by a primary iron nucleus differs from the shower formed by a primary proton. To 
investigate the difference in the case at hand requires a more detailed study.[] Note that the 
BH cross section plotted in Fig. 3 is inclusive over the mass of the BH. A significant portion 
of that cross section is due to the production of not so heavy BH's, through scattering with 
partons with small x. Moreover, as discussed above, the cross section a}jh is of the same 
order as the elastic gravitational scattering at ?/ ~ 1. In the latter processes a significant 
fraction of the neutrino energy is transferred to a few proton fragments. We then expect 
the resulting shower to resemble those induced by SM physics. In order to assess how well 
can one distinguish gravity induced showers from SM showers requires to take into account 
all these facts. This is an important point: if an excess of horizontal shower is observed, 
the shower shape information will be crucial to secure that the excess is not due to an 
underestimate of the (unknown) neutrino flux. 

Finally, in order to establish which process (gravitational elastic scattering, or black hole 
production) dominates the signal, one needs some knowledge about the energy dependence 



^For the observation of horizontal air-showers from neutrino primaries, see |27 
^See ra. 
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of the incoming neutrino flux. We have already estabhshed that BH production dominates 
energy loss. However, as the eikonalized cross section grows with E, if the neutrino flux J{E) 
decreases with E slowly enough, the number of BH events at energy E may be overshadowed 
by soft scattering events due to neutrini with energy ^ E. The signal is the number dN{E) 
of showers with energy between E and E + dE. In terms of the neutrino flux J{E) and the 
differential cross section we can write 

dE Je E' ^ ' dy^ E'' ^ ' 

Emax represents the energy at which (Jhh becomes larger than the inverse horizontal depth 
Xft"^. Neutrini with E > Emax interact right away and cannot generate horizontal showers. 
We have studied the above integrand by assuming J{E) ~ E". We found that, in the cases 
of interest, already for a < —2 the signal is dominated by events with large y, and then by 
black holes. (More precisely we flnd that the critical a's for n = 2, 3 and 6 are respectively 
equal to —1.76, —1.65 and —1.48.) This condition is satisfled for the cosmogenic neutrino 
flux in flg. 1 of ref. ||26| for which a ~ —3. If the cosmogenic neutrini dominate the flux, then 
black hole production and gravitational scattering at ?/ ~ 1, and not the softer processes, 
will dominate the signal in horizontal air showers. 

To conclude, we hope to have convincingly established that neither higher-dimensional 
graviton-mediated neutrino-nucleon scattering nor black hole production in TeV-gravity 
models can explain the observed cosmic ray showers above the GZK limit. Nevertheless, 
horizontal air showers may probe these scenarios. In this case, black hole production and 
gravitational deflection by a large angle will be the processes that dominate the signal. 
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